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We use the Kubo formalism to evaluate the contribution of acoustic 

O 
O 

phonon exchange to the frictional drag between nearby two-dimensional elec- 



^ . tron systems. In the case of free phonons, we find a divergent drag rate (r^J 

However, r^, 1 becomes finite when phonon scattering from either lattice im- 
perfections or electronic excitations is accounted for. In the case of GaAs 
quantum wells, we find that for a phonon mean free path smaller than 
a critical value, imperfection scattering dominates and the drag rate varies 
as ln(£ p h/d) over many orders of magnitude of the layer separation d. When 
^ p h exceeds the critical value, the drag rate is dominated by coupling through 
an electron-phonon collective mode localized in the vicinity of the electron 
layers. We argue that the coupled electron-phonon mode may be observable 
for realistic parameters. Our theory is in good agreement with experimental 
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results for the temperature, density, and (^-dependence of the drag rate. 
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I. INTRODUCTION 



Interactions between particles are a cornerstone of much of today's research in physics. 
In nuclear and high-energy physics, the effects of these interactions can be probed directly 
through scattering experiments. In condensed matter physics, inter-particle interaction ef- 
fects are enriched by the close proximity of other particles giving rise to a plethora of 
fascinating phenomena. However, direct measurement of these interactions in a condensed 
matter system is often a more difficult exercise, because of the indirect way in which scat- 
tering amplitudes are related to observables. 

Some time ago, Pogrebinskii and later PriceS proposed the following direct probe of 
interparticle interactions through a transport measurement. Place two 2-dimensional (2D) 
electron films close enough together and draw a current in one film. Through interlayer 
interactions, net momentum is transferred to electrons in the adjacent film, inducing a 
current there which can be measured. Due to technological difficulties in contacting the 
individual layers, decades passed before the first frictional drag experiment between 2D 
and 3-dimensional (3D) layers was performed.! The first experiments on this phenomenon 
between two 2D systems, as originally envisaged in Ref. [I] were performed by Gramila et 
al. for two electron layers,@S and by Sivan et al. for an electron-hole system.il In these 
experiments a current is drawn in the first layer, while the second layer is an open circuit. 
Instead of a current in the second layer there will be an induced electric field that opposes 
the "dragging force" from the first layer. The transresistivity p 2 i is defined as the ratio of 
the induced electric field in the second layer to the driving current density in the first, 

p 2 i ■ Ji = E 2 . (1) 

The stronger the interlayer interaction, the larger the magnitude of the transresistivity. (In 
this paper, we shall treat isotropic systems at zero magnetic field, hence P21 is diagonal). 
The transresistivity is often interpreted in terms of a drag rate which, in analogy with a 
Drude model, is defined by r^ 1 = p2\n\e 2 /m* where n\ is the electron density of the driving 
layer and m* is the electron effective mass. 
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These experiments spurred a large body of theoretical work both on electron-hole 
systems! and on electron-electron systems. EH! Most of this work has focused on interlayer 
Coulomb interaction, the most obvious coupling mechanism and the one considered in the 
original theoretical papers.0 However, it was clear from the start that the experimental re- 
sults were inconsistent with a purely Coulomb interlayer interaction, which predicted a low 
temperatureiBEl {ksT <C ef,i, £f,2, where e Fj i is the Fermi energy for layer i) transresistivity 
of the form 

pn = (_ M C(3)vr {k B Tf 1 1 (2) 

21 V e V 32 e F<1 e Fj2 (k F)1 d) (k F;2 d) {q TF d) (g T F^) ' 

where ( is the Riemann zeta function, d is the interlayer separation, kp,i is the Fermi wavevec- 
tor for layer i, and q^F is the Thomas Fermi screening wavevector of the 2D electron gas. 
This expression is based on the random phase approximation (RPA) for the screened inter- 
layer Coulomb interaction and applies for q^pd 3> 1 and k F ^d 3> 1. From inspection of Eq. 

(0) , one notes three important characteristics of the Coulomb drag for low temperatures: 

(1) the scaled transresistivity p 2 i(T)/T 2 is a constant; (2) p 2 i is a monotonically decreas- 
ing function of the density of either layer (so long as k F ^d ^> 1); and (3) p 2 \ oc d~ 4 . The 
experimental results at around 2 — 3 K, on the other hand, showedBSEl (1) a well-defined 
peak in the p 2 i(T)/T 2 , (2) a local maximum in p 2 \ near equal layer densities and (3) an 
approximately d- independent residual p 2 \, after subtraction of the Coulomb contribution 
which can be identified experimentally by its simple d and T dependencies. Furthermore, 
the experimentally measured magnitude of p 2 i was generally larger by about a factor of 2 
than the value predicted by the Coulomb interaction alone. Another momentum transfer 
mechanism was clearly involved. 

From the outset it was understood that exchange of acoustic phonons was the most likely 
candidate for this second momentum transfer mechanism. Exchange of phonons often dom- 
inates electron-electron scattering contributions to the resistance of bulk metals.^ The peak 
in the temperature dependence of p 2 i(T)/T 2 is reminiscent of features in the temperature 
dependence of the acoustic phonon limited mobility, and occurs around the Bloch-Gruneisen 



temperature Tbg = 2k^ 1 hcikp associated with the acoustic phonon modes, (q is the longi- 
tudinal acoustic phonon velocity.) The phase-space available for scattering is largest for 2k p 
transfers, partially explaining the enhanced drag when kp,i = &f,2- Finally, the long-ranged 
exchange of phonons between electrically isolated systems is not an unknown phenomenon. 
Exchange of phonons between two 3D systems separated by ~ 100 /im has been observed 
previously^ and related effects are expected to be observable in a superlattice if the driving 
electron layer is The theoretical challenge is to explain the magnitude of the observed 

drag and its dependence on layer separation and density. 

Since interactions of acoustic phonons with electrons are relatively weak in GaAs, the fact 
that phonon-mediated and Coulomb contribution to the drag are often comparable seems 
mysterious. However, we show below that the obvious calculation, in which a free phonon 
propagator substitutes for the Coulomb interaction, leads to a divergent drag resistivity. 
The large but finite drag rates which are observed experimentally can be explained in terms 
of scattering and interaction effects which alter the phonon propagator. 

Despite the apparent importance of phonon exchange in drag measurements, there has 
been less theoretical work on this mechanism than on the Coulomb coupling mechanism. An 
incoherent phonon exchange model studied by Gramila et al. produced a drag rate which was 
too weak to account for the observed transresistivity.S Tso et aM addressed the question of 
whether exchange of virtual phonons could give a contribution strong enough to explain the 
magnitude of the transresistivity, but did not use an electron-phonon coupling model which is 
realistic for GaAs/Al x Gai_ x As systems. Zhang and Takahashi0 included all relevant phonon 
contributions; however, due to an incorrect effective two-dimensional electron-phonon inter- 
action, they predicted a short range phonon-mediated drag. In spite of the use of radically 
different models for the phonon-mediated process, both these calculations yielded a temper- 
ature dependence for the drag in reasonable agreement with experiment. The similarity is 
not surprising because of the common appearance of the Bloch-Griineisen temperature scale 
associated with the acoustic phonon mode. Comparable features in the drag, associated 
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with plasmon modes of the electronic system, were predicted00 and observedlM! at higher 
temperatures. 

In this paper we report on a detailed examination of the phonon exchange mechanism for 
drag, using a model which we believe to be quantitatively reliable for GaAs/AlAs quantum 
well systems. We address distance-, temperature-, and density dependence of the transre- 
sistivity. We show that two different regimes of layer separation dependence can occur, 
depending on the phonon mean free path and the electron-phonon coupling constant. We 
demonstrate the existence of a coupled electron-phonon mode, which for long phonon mean 
free paths leads to a large enhancement of the drag. 

The outline of the paper is as follows. In Section |I| we use a Kubo linear response 
formalism to obtain a formula for the transresistivity which is sufficiently general to permit 
the incorporation of a finite phonon mean-free-path and renormalization of the phonon 



propagator due to coupling to the electronic layers. In Section |T| we discuss the phonon 
mediated e-e interaction and explain how the relatively weak electron-phonon interaction 
can lead to a surprisingly large contribution to the drag. The two different regimes of £ p ^ 
are discussed in Section [TV] and |V|, where approximate expressions for p 2 \ are derived and 
special attention is paid to the layer separation dependence of the drag. Detailed numerical 



results are presented in Section [VI] before we conclude in Section |VII| with a summary of 
results. 



II. FORMALISM 

Identical theoretical expressions for the drag rate due to Coulomb interactions been 
obtained in several different ways. The most physically transparent derivation is based on 
semiclassical transport theory.0 More elaborate fully quantum mechanical derivations based 
on memory function' or Kubo formula approaches" yield identical results for the large 
kpi limit at zero magnetic field but are more flexible and, in particular, can be applied in 
the presence of an external magnetic field. (Here I is the electronic mean free path.) We 
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show below that at lowest nonvanishing order in the electron-phonon interaction, phonon 
exchange yields an infinite result for the drag rate. The Kubo formula approach, which 
we use in this paper, is most convenient when interlayer interactions need to be treated 
beyond leading order, because of the powerful diagrammatic perturbation theory expansion 
available to evaluate the influence of interaction terms on the appropriate current-current 



correlation function. The required calculation is an adaptation of those described in Refs. |13 



and 15 and is outlined below. 



A. Electron-phonon interaction Hamiltonian 

We consider frictional drag between two GaAs quantum wells. We define the plane of 
the quantum wells as the x-y plane. The distance between the centers of the two quantum 
wells is d, and the width of the two wells is L. We assume that the electron number densities 
in each well, n\ and ri2, are such that only one subband of the quantum well is occupied. 
The formalism can easily be generalized to accommodate more occupied subbands. 

We include interlayer and intralayer Coulombic electron-electron interactions and the 
coupling of electrons in either layer to the three-dimensional (3D) phonons of the semicon- 
ductor host.EI We denote 3D wavevectors by upper case letters and their projection onto 
the x-y plane by the corresponding lower case letter so that Q = (q, Q z ). The electron 
creation (annihilation) operator in layer i is cf(k) (cj(k)) with implicit spin indices, the 
phonon creation (annihilation) operator for polarization A is d\ q (oa,q)j and the subband 
wavefunction of electrons in well % is tpi(z). With these definitions, the electron-phonon 
interaction contribution to the Hamiltonian is given by 

tfe-ph = V- 1 ' 2 E E E M A(q, Q g ) A^q. pi(-q)Fi(Q z ) (3) 

A i=l,2q,Q z 

where V is the normalization volume, 

/oo 
dz \if t (z)\ 2 e-^ z (4) 

A\,Q = Oa,Q + ^A,-Q 
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ft (q)=E4(k)Q(k + q) (5) 

k 

and M\(Q) is the bulk electron-phonon coupling constant. At temperatures much lower 
than the Debye temperature one can neglect the Umklapp process in the electron-phonon 
interaction Hamiltonian. 



B. Kubo formula transconductivity 

The Kubo formula for linear response offers an expression for the transconductivity tensor 
which is defined by 

J 2 = £ 21 Ex. (6) 

In the absence of magnetic fields, J 2 and E x will be anti-parallel, and a 2 i is a diagonal 2x2 
tensor. 

The derivation sketched below for a 2 i is very similar to the one given previously in 



Ref. 15, in which the reader can find further details. The transconductivity is given by 



o-%(k,n) = —n a 2 r e \k,n) (7) 

where n 2 7' re *(k, £7) is the Fourier transform of the retarded current- current correlation func- 
tion, and a and 7 are Cartesian indices. The retarded correlation function is evaluated by the 
standard analytic continuation of the (bosonic Matsubara frequency) Fourier components of 
its imaginary time counterpart .El The imaginary time correlation function is calculated in 
perturbation theory 

11 21 (x - x , r - r ) - , (8a) 

(8b) 

where (■ ■ -)o denotes a non-interacting system thermal average, and T T {- ■ ■} is the usual 
r-ordering operator. The S'-matrix is expanded in powers of the interaction Hamiltonian 
and the resulting non-interacting system correlation functions are evaluated with the aid of 
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Wick's theorem. As usual, the denominator in this expression cancels the "disconnected" 
terms in the diagrammatic expansion. 

If we include only electron-phonon interactions for the moment, the lowest non- vanishing 
term appears at fourth order and gives the following contribution to the current-current 
correlation function: 

U-2i (k = 0,i£l n y A > = dr dr\ \ dr 2 / rfr 3 / dr 4 exp(zfi n r) 

4 Jo J o J o Jo Jo 

E E E E M Al (Q 1 )M A2 (Q 2 )M A3 (Q 3 )M A4 (Q 4 ) 

Ql,Al Q2,A2 Q3,^3 Q4.A4 

(T T J?(q = 0, r) pi(-qi, n) pi(-q 2 , r 2 ))o 
(^r jf (q = 0, 0) p 2 (-q 3 , r 3 ) p 2 (-q 4 , r 4 )) 

( T r ^Q 1 ,A 1 (n)^Q 2 ,A 2 (r2)^Q3,A 3 (r3)^Q 4 ,A 4 (T-4))o (9) 

where A is the 2D system area. 

The Wick's theorem factorization of the phonon operator product expectation value leads 
to the product of two bare phonon Green's functions, defined by 

7jf(Q, t-t>) = -(T r {i AiQ (r)i Aj _ Q (r')})o. (10) 
It follows that, to leading order in electron-phonon interactions, 

n 2i 7 ( k = o, zfi) = — \o E E A 2 (-q, -q, ~ iuJ - -«*>)A?(q, q, «*> + ^, «*>) 

X / -% > /M(A)/->(-(;OIJ/ A (q.(,MriA (q.(A./^) 



/ f| E ^i(^)^(-Q z )|M A (q, Q 2 ) | 2 Df (q, Q 

A 

2nh 4; 



/ f§ E Fi(g;)F 2 (-g;)|M A ,( q , Q'j| 2 £>S?(q, Q' z , <n + «*) (n) 



where 

A(q, q; 2W n , iwj = -A / dr ± / dr 2 (T T j(q = 0, 0)p(q, n)p(-q, -r 2 )> 

Jo Jo 

exp(i^ n Ti) exp(zu^T 2 ) (12) 

The Feynman diagram corresponding to this contribution to the correlation function (|TT|) is 
shown in Figure |l|. 



It will turn out to be important to account for disorder and anharmonicity in the lattice 
system. We will do so using a phenomenological approach by introducing a phonon mean 
free path, i^. However, since the intrinsic bulk phonon mean free path may exceed the 
dimensions of the sample, we should in principle take the surface scattering explicitly into 
account. For simplicity we will nevertheless use a single phenomenological mean free path 
and later, when we discuss the long mean free path limit in more detail, take boundary 
effects into account. The phonon Green's function is therll 

D x (Q,iu n ) = -- 2UA,Q — ^— — — (13) 

(u n + {c x /2£ ph ) sgn(u n )y + uj z XjQ 



where U\,q = C\yQ 2 + Q% and l p h is the phonon mean free path. 

Note that the wavevector arguments in Aj are 2D, since Aj is a property of the 2D 
electron systems. One can therefore sum over A and integrate over Q z to obtain a phonon- 
mediated effective interaction 

Dy(q,iu; n ) = J H£ |M A (Q)| 2 F l (g 2 )F J (-g,) J D A (Q,^ n )- (14) 

A 

This effective interaction is the 2D Fourier transform of the product of the phonon propagator 
between the layers and the electron-phonon interaction in each layer. With this definition 
Eq. ( pTTj) becomes 

— 1 1 

n^i 7 (k = o, in) = 51 — J] ^2 (-q, -q; -iu - ito, -iuj)Aj(q, q; iu + iQ, iu) 
z An q rip iuJ 

T>2i(q,iuj n )'D2i(q,ift n + iuj n ) (15) 

This expression for n 2 i is the same as in Ref. |15|, except that the interlayer Coulomb in- 
teraction is replaced by the phonon-mediated effective interaction. From this point on, 
the formal steps are identical to the Coulomb case. Performing the summation over iu, 
continuing to real frequencies and taking the f2 — > limit, we obtainlil 



q 

>7 



dn B (u) 



du 

x A^— q, — q, — u — i8, —uo + i<5)A"(q, q, u + i5, u — i5). (16) 
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C. From Transconductivity to Transresistivity 



For further progress it is necessary to make some assumptions about the electronic sys- 
tems. We will assume that the 2D electron layers are good metals with large kpl where 
£ is the electronic mean free path. It can then be shown@0Jli that whenever the trans- 
port scattering time, r tr , is independent of energy, the function A is related to the electron 
polarization function x(q, to) by 

m* 

Here m* is the electron effective mass. The relation (|l7]) is a property of the 2D electron 
layers only, and is not dependent on the phonon degrees of freedom. Assuming further that 
| c"i2 1 <C On the transresistivity can be approximated as follows: 

P21 = ~ • (18) 

°"ll cr 22 ~~ cr 12C r 21 C r ll cr 22 

Using the Drude expression (an = e 2 niT tr Jm* ) for the intralayer conductivities, which is 
valid under the above assumptions, the transresistivity due to the electron-phonon interac- 
tion is given by the following explicit expression, 

P21 = ~r- 2 T^-TZ^ / T~ ^21 (q,^) . 7J7Z /0 , — • 19 

Ae l nin2k B T A ^ J-oo 2ir smh (nuj/2k B T) 

The same expression for the phonon-exchange contribution to the drag can also be derived 
using semiclassical Boltzmann transport theory and a collision term with transition matrix 
elements calculated by summing over virtual and real intermediate states with absorbed 
and emitted phonons. Ambiguities can arise in that approach, however, from the portion of 
phase space where the relevant energy denominators approach zero. As we explain below this 
part of the phase space is important in determining the drag resistivity. Our Kubo function 
derivation allows finite phonon mean-free-paths, which remove any spurious singularities, to 
be incorporated into the calculation in a consistent and unambiguous manner. 

Note that the explicit dependence of the transconductivity on the transport lifetime is 
absent in the transresistivity. (An implicit dependence remains through the dependence 

11 



of the polarization function on disorder.) This aspect of the final expression has been 
emphasized in recent work of Swierkowski, Szymahski and Gortel,0 is not accidental, and 
emerges naturally in a force-balance approximation where the drag force simply cancels the 
rate of momentum transfer per particle from the current carrying layer to the open layer. 



D. Coulomb interaction and screening 

In reality, electron-electron and electron-phonon interactions are simultaneously present 
and both should be included in a drag calculation. To leading order, Coulomb interactions 
can be incorporated by simply addingi! the interlayer Coulomb interaction to the phonon- 
mediated effective interaction in Eq. (|I9|) . One class of higher order terms in which intralayer 
Coulomb interactions appear is captured by replacing the electronic polarization functions 
which appear in Eq. (|19|) by their interacting system counterparts. However, the most 
essential higher order terms are those which account for the screening of both phonon- 
mediated and Coulombic interactions. In the RPA, the total interlayer screened interaction 
is given by 

W^ M = V ^ + . U ^\ (20) 
where Uij(q) is the unscreened Coulomb interaction and 

e(q, u) = [l- (p n + U u )Xi] [1 " (^22 + U 22 ) X2 ] - (Z> 21 + U 21 ) 2 XlX 2- (21) 

(See for example Ref. 0. The form of the phonon-mediated interlayer interaction in this 
reference is incorrect, however.) e(q,u) is the effective dielectric function for interlayer 
interactions in the RPA. Notice that when the contribution of interlayer interactions to the 
screening can be neglected, e is simply the product of the dielectric functions for the two 
layers, corresponding to independently screened electron-phonon interactions in each layer. 

Coulomb and phonon-mediated interactions can be simultaneously included in the tran- 
sresistivity simply by replacing V by V(/ total in Eq. (|T9j). Note that the transresistivity is 
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not strictly the sum of purely "Coulomb" and "phonon" contributions, since there are in- 
terference terms proportional to U x T> in the |W|° tal | 2 . However, the Coulomb contribution 
is large only for q < 0.5kp, {U(q) = 27re 2 exp(— qd)/q when the finite thickness of the elec- 
tron layers is neglected, whereas contributions from the D term come predominantly from 
q pa 2kp. Hence, the interference terms are usually negligible, and in practice we will treat 
the Coulomb and phonon contributions as if they were incoherent. In what follows, we will 
concentrate on the "phonon contribution," which will be calculated from Eq. (|l^), with only 
the T> term in the numerator, 

W 2 i(g,w) = ^i (22) 

Note that it is important to retain the coupling between the Coulomb and phonon terms in 
e, since this can influence the q ~ 2k p contribution to the transresistivity. 

III. THE PHONON MEDIATED INTERACTION IN GAAS/AL x GAi„ x AS 

SYSTEMS 

As mentioned earlier, the acoustic electron-phonon interaction is weak in GaAs and 
ALGa^^As. In this section we discuss quantitatively what "weak" means, and explain how 
phonons can give an important contribution despite this weakness. In the following section 
we discuss how and when screening affects the phonon mediated interaction. 



A. Electron-phonon coupling in GaAs/AlGaAs 

In GaAs/AlGaAs systems, electrons couple to acoustic phonons via deformation po- 
tential and piezoelectric couplings. Since we are concerned with low energy excitations, 
only acoustic phonons in the long wavelength limit have to be considered. In this limit, 
the squares of the e-ph coupling strengths for longitudinal and transverse phonons are,0 
respectively, 
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|M(Q)| 2 = ^ 

2 get 



(23) 



where g is the mass density of the crystal, D is the deformation potential, ehu is the 
piezoelectric constant, the c\ are sound velocities for longitudinal and transverse phonons, 
and the A\ are the anisotropy factors,^ 

MQ) = ^ (25a) 



,2 



B. Approximate analytic form of T>(q, u) 

While it is possible to obtain exact expressions for T> from Eq. (|14D including the full 
anisotropy functions (Eqs. (|25a|) and ( |25b|) ) and the form factors for infinite square wells, 
these are extremely complicated. Therefore, we shall make some well-controlled approxi- 
mations explained below which do not significantly affect the final results for the computed 
transresistivity. 

As we shall see shortly, the phonon-mediated effective interactions are important only 
when u} is close to c\q. For these u's the integral over Q z in Eq. (TfiJ) is dominated by 
contributions from near Q z = 0. We therefore remove all Q z dependent factors, except for 
the small energy denominator, from the integral which defines T>. For example, one can set 
Q z = in the anisotropy factors given in Eq. ( |25a| ) and ( |25b| ), yielding 

At « 0, At « -. (26) 

Then the phonon-mediated effective interaction for two equivalent infinite square wells with 
width L and center-to-center separation of d is 



vq Vij(q, u) « -^rrf-tij - C DP Bij ( qdy/1 - Cf\ qLyJ 1 - Q 

-C PE , kF By (Ww^ \ I cr ) • (27) 
14 



Here, u = m* /ttTi 2 , the 2-dimensional electron gas density of states (so that uqD is dimen- 
sionless) , 

b = ^ + A~ (») 



c x q 2q£ ph 
(eh u ) 2 m* 

D 2 m*k F 
2nh 2 c 2 g 



C PE = (29) 
Cop = ^ (30) 



and 



Bii(x,y) = { y2+ * 2 W y fy^ 2[e L) > ' l ~ h (31) 
exp(-x) (^) 

In the expressions above, the square root with a positive real part should be taken. 

Inserting numerical values for GaAs (m* = 0.067m e , q = 5.14 x 10 5 cm/s, q = 3.04 x 
10 5 cm/s, q = 5.3g/cm 3 , ehu = 1.2 x 10 7 eV/cm and D = — 13.0 eV), gives the following 
dimensionless coupling constants: 

1 6 cm -1 

C PE w 1.64 x 10- 3 x (32) 

C DP « 2.7 x 10- 3 x _ 1 (33) 

10° cm 1 

For the same subband wavefunctions the bare interlayer and intralayer Coulomb inter- 
actions are given by 

Uij(q) = — B i:j (qd, qL), (34) 

where q^F — 27re 2 t / is the Thomas Fermi wavevector. (We have absorbed the bulk dielectric 
constant of the semiconductor in the electron charge). 



C. Strength of £> 2 i 

For the RPA screened Coulomb interaction (C/rpa) in a single layer system, the magnitude 
of z^o^rpa approaches 1 at long wavelengths. Although the corresponding value for double- 
layer systems is smaller,! it is useful to compare v§T>2\ to this value. For our present 
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illustrative purpose we concentrate on the deformation potential term (i.e. the D 2 term 
in Eq. (|23|)), which turns out to dominate except for very low density electronic systems. 
Then, in the limit of vanishing quantum well widths and £ p h — > oo, we obtain 

v U 21 (q,u + irf) « -C DP = = exp ( -d\jq 2 - u 2 cj 2 J (35) 

qk F c 2 ^l-0J 2 /{qci) 2 v ' 

Note that the magnitude of the effective interaction diverges as q — > u/ci from above or 
below. (We point out that the transresistivity is strongly dependent on the width L, and 
therefore one should not take the L — > limit when computing T>. We use the full form 
given by Eq. (p7|) in subsequent numerical calculations.) 

The small prefactor C^p guarantees that the phonon- mediated effective interaction is 
small compared to the Coulomb interaction except near uj ~ qc\. The large value of the 
interaction in this region of phase space reflects the large phase space for intermediate states 
with small or vanishing energy denominators when u ~ qc\. (Note that the phonon energy 
varies slowly with Q z for Q z near zero.) 

The importance for drag of the sharp peak in D 21 (q, uj) is enhanced by the fact that it 
appears squared in Eq. (|T^). As a result, phonons do play an important role in frictional 
drag, even though the typical value of z/ X> 2 i(q, w) is small. In fact, if we ignore the effects 
of screening and let £ p h — > oo, we get an infinite transresistivity. This is easily seen from 
Eq. fl35|) . The absolute value £> 2 i diverges like (c\q — uj)^ 1 ^ 2 as u approaches c\q from 
both above and below. Inserting this bare form of the interaction into the Eq. ( |I9|) gives 
a I c a9 — divergence in the energy transfer integral for every q, yielding an infinite 
transresistivity. This point does not seem to have been emphasized in the existing theoretical 
literature on this subject. As we show in the next two subsections, including either a finite 
mean free path or screening the interaction dynamically removes this spurious divergence. 

In Fig. ^ we plot h / \V 2 i(2k F , u)/e(2kF, u))\ near the longitudinal resonance for different 
values of £ p h- It is useful to compare ^ol^iAl to the screened interlayer Coulomb interaction 
^o^21,rpa oc exp(— qd)/(l + <?tf/<?)- In contrast to the phonon-mediated interaction near the 
resonance, ^oL^i.rpa decreases rapidly below unity as a function of well separation. Fig. |2|, 

16 



partly explains how the "weak" phonon-mediated interaction can compete with the Coulomb 
interaction as a mechanism for drag. 



D. Effect of Screening 

From Eq. (|27|), one sees that the presence of a finite £ p h cuts off the divergence and leads 
to a finite transresistivity. Since the divergence in the integrand is of the form l/\u — c\q\, 
it follows that p 21 would be proportional to ln(^ p h) if screening were not important. When 
screening is included the drag resistivity does not diverge. 

As shown in Eqs. ( |21~D and (fZUD , screening is accounted for by dividing the bare interlayer 
interaction £>2i(<?, ^) by a dielectric function e(q, to). Whenever T>2i(q, uj) diverges at uj = c\q, 
so does e(q,u). The screened interaction is therefore nondivergent, even when — > oo. 

As we explain below, screening becomes important for the phonon-exchange drag only if 
the phonon mean free path exceeds a critical value £ p h, C mt- Since £ p h, C rit is close to realistic 
values, we investigate the two regimes separately in the following two sections. 



IV. LONG MEAN FREE PATH LIMIT 

In this section we focus on the large phonon mean free path limit, where a coupled 
electron-phonon mode turns out to be of utmost importance. First we discuss the ideal case 
of infinite £ p h, where the coupled mode is broadened by the coupling to the electronic system 
only. In this case we find an analytic form for the coupled mode contribution. Secondly, in 
Section [IV B| , we discuss how the collective mode contribution is modified by a finite mean 



free path and the conditions for its experimental observation. 



A. Infinite phonon mean free path limit 

The approximate analytic results discussed below apply only for layer separations smaller 
than a large but finite maximum value which we specify below; for still larger layer separa- 
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tions the analytic analysis is less revealing and we have relied more strongly on numerical 
studies. 

We show below that for £ p h — ► oo, the real part of e(q,u) vanishes and the imaginary 
part is small enough to yield a sharp collective mode for uj just below qci which contributes 
strongly to the frictional drag between the two layers. To make the following discussion 
as transparent as possible we limit our attention to the case of identical electron layers so 
that V 2 2 = T>n and U\\ = £^22 • We locate the collective mode frequency u>o(q) by solving 
the equation Re[e(g, ujq)] = 0. We will be interested in momentum transfers near 2kp which 
give the main contribution to the drag and systems with kpd ^> 1 so that we can neglect 
U2i{q) oc exp(— 2kpd). The small layer separation approximation mentioned above consists 
of setting V u m £> 12 = V in Eq. which is justified for 

d < Tar I — = dB ^ 
as we show later. With these assumptions 

W 2lM « - - - (37) 

(1 - Unx){l - Unx - 2£>x) 

The collective mode in which we are interested occurs near uo = c\q vpq and hence is in 
the low frequency regime for x{q, w)- Here vf is the Fermi velocity of the electronic system. 
We therefore approximate Rex(<7, oj) ~ xil^ = 0) = — u . We write the imaginary part of 
the polarization function in the form lmx{q, uj) = —VqS. In the low frequency limit of the 



RPA, 5(q, ciq) ~ q/ (vpyjl — q 2 /4:k F ), except near 2kp where it is approximately Jci/2vf. 5 
is dimensionless and has a wavevector and frequency dependence which will be left implicit 
in the following discussion except where emphasis is important. However, it is important 
for the following discussion that 5(q, ciq) is generally small compared to one. 

Because we are interested in large wavevectors near q = 2kp we include only longitudinal 
deformation potential coupling. The calculation including transverse piezoelectric mode 
follows mutatis mutandis. From the previous section, we have, for uo close to but smaller 
than aq, 



voW2i(q,u) = — , =. (38) 

e (!fejl - uy^q* - 2C DP ) - i2C DP 5 



In this equation eo = 1 + <?tf/<7 is the RPA static dielectric function for an isolated layer. 
The collective mode occurs where the real part of denominator of Eq. (|38|) vanishes at 
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\ V k F e 

This collective mode results from coupling of the electron layers to phonons with u close to 
ciq. Expanding the denominator of Eq. ([38]) around the pole we find that for u close to ujq 

Hw 2l {qM\ 2 = c 4^fi ( 4 °) 

e*k 2 F 5 [0J - ^o) 2 + 

where the width of the collective mode resonance is given by 

r 4~5 Cl q 3 C 2 DP 



,3p 



(41) 



Note that the width of the resonance is small compared to the shift of the resonance from 
Qq and that the resonance line shape is approximately Lorentzian only if <5 <§C 1 . Where this 
condition is not satisfied Eq. ([40]) will not be accurate. 

Numerical calculations discussed later demonstrate that for £ p h — > oo the drag is domi- 
nated by coupling associated with this collective mode resonance. Since the effective inter- 
action has a more rapid frequency variation than other quantities in the expression for the 
transresistivity (Eq. flT9D) we may approximate the screened interaction near the resonance 
by a 6 function: 

WoW 21 (q,u)\ 2 = ^^^(u-cuo) (42) 

This approximation allows the frequency integral in Eq. (|T9|) to be performed and the con- 
tribution from the coupled mode to the transresistivity may be expressed as 

fh\ 7iK x f°° , (q a \ Jy lmx{q,c x q) 
p21 ^[7i)nWJ dq {4) smh\hc x q/2k B TY (43) 

where a = 6 for the longitudinal phonons and a = 2 for the transverse phonons, and where 
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(Air) 3 cfg 2 
h(eh u ) 4 

Kt = 2*%W ' (44b) 

(44c) 

The q 6 /el term in the integrand above implies that for T <^ Tbg = 2hcikp/kB, the main 
contribution to the integral comes from the q = 2kp region; i.e. large-angle scattering 
dominates the phonon mediated drag. 

Comparing with Eq. (|2]) we see that, at least for £ p h — > oo, the phonon mediated drag 
can be comparable to or stronger than Coulomb drag, in spite of the weak electron phonon 
interactions. Crudely the condition which needs to be satisfied is that Tbg is low enough 



or the layer separation is large enough that CppyJ Ci/vp ~ {k B T-Q G / '€p) 2 (qTFd)~ 2 (kFd)~ 2 . 
For typical layer densities in GaAs this condition is satisfied for layer separations larger 
than a few tens of nanometers, consistent with experimental observations. When phonon 
mediated drag is dominant, P21 will be proportional to temperature for T ^> Tbg- Below 
Tbg there will be a crossover to a regime where the piezoelectric contribution dominates and 
the temperature dependence goes approximately as T 5 . This can be seen from Eq. (^) by 
defining a new integration variable proportional to q/T and taking the g-dependence of eo 
into account. We shall discuss this further in Section |V1 A| . At extremely low temperatures, 
the assumption that the drag rate is determined by the resonance u ~ c\q will break down. 
P21 is then dominated by the u = limit of T> 2 ± and will revert to the familiar T 2 law for 
carrier-carrier scattering in a Fermi liquid. 

Recall that we used T>2i ~ to obtain the above results. From Eq. 

i> !i(?> w o ) j(9)) 



m^M)'^ { ~ 2qd ^ 2{1 ~^ My (45) 

Hence, for the approximation V^—DW to be valid, one must have 2qdJ2(l — cu /ciq) 



1. Together with Eq. (|55D, this gives the condition Eq. ([36]). For £ p h — > oo, the preceding ap- 
proximate calculation implies that there is no layer separation dependence until this length, 
which for GaAs and typical densities corresponds to d ~ 5000 A, is reached. Numerical 
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results show, however, that there is a weak distance dependence on the transresistivity in 
this large I regime. This stems from the presence of a relatively long non-Lorenzian tail for 
uj < Qq which contributes significantly to the integral, as we will discuss below. 

At larger layer separations, the interlayer phonon propagator V 2 i is reduced at the col- 
lective mode frequency and the two electron layers interact with the phonon system more 
independently. We must then use a more refined expression for the effective interaction: 

^" (i-K.W-iftx- (46) 

A competition occurs between a decline in the coupling due to suppression of the interlayer 
propagator D 2 i and resonant enhancement of the electron phonon interaction near each 
electron layer. 



B. Coupled electron-phonon mode with finite mean free path 

The above analysis is based on phonons with an infinite mean free path, apart from 
the finite lifetimes due to interactions with the electronic layers. Any real system will have 
imperfections which will make the phonon life time finite even when the electron layers are 
not present. Even when the lattice is perfect and free of isotopic impurities, anharmonicity 
and boundary scattering will cause phonon modes to decay. In the following we represent 
all these effects in the simplest possible way by assigning a common phenomenological mean 
free path £ ph to all modes. 

The occurrence of coupled phonon-plasmon collective modes, signaled by by a zero in 
the Re[e(g,u;)], requires a cancellation between phonon-mediated and Coulomb interaction 
contributions. In the present section we will consider the limit in which the phonon-mediated 
contributions to e(q,u>) can be neglected, i.e. the limit in which \vqD\ <C 1 for almost all 
energy and momentum transfers and for both intra-layer and inter-layer propagation. We 
will see that this condition is satisfied except at long phonon mean free paths. 

For uj close to qc\ and finite £ p h (but large compared to q^ 1 ), 
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|j> z>2i| « (qC DP /k F ) l^i-Cfr 1 - (47) 

where C/ is defined in Eq. (p8|), and 1 — Q ~ 1 — (u/qci) 2 + i/q£ p h- There is a critical 
value £ p h,crit such that 1 — (co> /qg) 2 < l/<Z^ p h,crit> in which case | Refill never becomes large 
enough for Re[e] = to have a solution and the collective mode ceases to exist. This critical 
value is given by 

egfcj. (1 + ?tf/2M 2 
Writ- 4C 2 p? 3 32C£ P A* ' 1 j 

where, since large momentum transfers dominate the drag, we have set q = 2k p. £ p h, C rit 
is even longer at smaller g's and it seems improbable that the collective mode discussed 
in the previous section would ever be evident in inelastic light scattering studies of long- 
wavelength electronic excitations. In GaAs, and for kp ~ 10 6 cm _1 , the critical mean free 
path is approximately ^ p h,crit ~ 0.2 mm. The actual mean free path of course depends on 
the sample in question. 

It is often the case that scattering of phonons off the boundaries of the sample can be 
accounted for by taking the mean free path in the absence of bulk scatterers to be equal 
to the sample size. We argue below that for the present problem, the sample size in the z- 
direction is usually irrelevant and that the maximum mean free path is given by the typically 
larger lateral dimension of the sample. 

The phonons which contribute to the coupled mode are those which are confined around 
the two-dimensional electron gases. The extent of the confinement is given by range T>2i in 
the z-direction {i.e., perpendicular to the 2DEGs). To study the ^-direction range of the 
phonon field participating in the collective mode, we look more closely at which z-region 
actually contributes to the coupling. We write the effective phonon mediated interaction 
given in Eq. fll4|) as 

D 2 i{q,w + i6) = J dzdz' \ip2{z)\ 2 \ipi{z')\ 2 K{q,uJ + i6, z - z'), (49) 

where K(q, u + iS, z — z') is the Fourier transform of the phonon Green's function and the 
electron-phonon coupling matrix element with respect to Q z , 
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K(q,cu + iS,z)=J2 J ^-e tQ * z D x (Q,LU + t5)\M x (Q)\ 2 . (50) 

A 

For a given q and u, the spatial extent in the z-direction of K(q, uo + i8, z) gives the range 
of the phonon field. Since the collective mode frequencies are very close to uo = Qq, the 
relevant phonon wavevectors in the z-direction are small, and the same approximations used 
previously can also be applied here, yielding an approximate form for K given by Eq. (|35|) 
with d replaced by \z — z'\\ i.e., K(q, u + 5, z — z') ~ exp(— \z — z'\Jq 2 — uj 2 /c}). The spatial 
extent of the phonon field participating in the collective mode can be found by substituting 
the frequency for mode, uq, into this exponential form, giving K ~ exp(—\z — z' \ / ds) , where 
ds was defined previously in Eq. (j36|). For GaAs at n = 1.5 x 10 11 cm -2 ds is approximately 
5000 A. Any boundaries or imperfections beyond this range in the z-direction have negligible 
effect. 

It is plausible that some existing experimental results are for samples with £ p h in the 
lateral direction comparable to this critical mean free path and hence partially reflect col- 
lective excitations of the electron-phonon system. However, more experiments are needed 
before definite conclusions can be drawn. 



V. SHORT PHONON MEAN FREE PATH LIMIT 

In the present section, we will consider the limit in which phonon mediated contributions 
to e(g, uj) can be neglected, i.e., £ p h £ p h,crit- We begin our discussion by addressing the 
distinction between the real and virtual phonon contribution. This we do partly due to 
somewhat confusing use of these terms in the existing literature. 



A. Real and virtual phonons 

The phonon mediated effective electron-electron interaction D2i(q,uj), obtained here 
from diagrammatic perturbation theory, can also be derived (with a little more work!) using 
elementary time-dependent perturbation theory starting from an expression of the formEl 
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Wf , u) . E MM), (51) 

j &i — &I + IT] 

Here \f), and \I) are the initial state, the final state in which momentum has been 
transferred between the layers, and an intermediate states in which the momentum to be 
transferred is carried by a non-equilibrium phonon. The infinitesimal imaginary part, it], 
in the denominator enforces causalityH and plays a crucial role when intermediate and 
initial or final states are close in energy. In time- dependent perturbation theory, this form 
of effective interaction determines the transition rate between initial and final states when 
the perturbing term in the Hamiltonian, the electron-phonon interaction in the present 
case, does not directly couple initial and final states. We therefore expect the phonon 
mediated effective interaction to play the same role in transport experiments as the interlayer 
Coulombic interaction which does have direct matrix elements between initial and final 
states. The intermediate state need not conserve energy and, in the jargon of time-dependent 
perturbation theory, is consequently referred to as a virtual state. 

The intermediate states are ones where a phonon of momentum Q = (q, Q z ) has been 
created or destroyed by an electron in layer 1 or 2. Summing over four intermediate states 
for each Q z and performing thermal averages gives the following expression for transition 
matrix elements with 2D wavevector transfer q 



V 21 (q,cu) = l d §^\M l {Q)\*F 1 {Q z )F 2 {-Q z 



1 1 



(52) 



U - UJq + 17] UJ + UJq- 17] _ 

Eq. (|TJj) reduces to this form for l p \ l — > oo except that we have, for simplicity, retained only 
the longitudinal phonons in the present discussion. 

Recalling that Two is the energy transferred between the layers and Hujq is the energy of 
the intermediate phonon, when the denominator vanishes in Eq. fl52|) energy is conserved 
in the intermediate state. Therefore, the real (imaginary) part of the term in the square 
parenthesis gives the virtual (real) phonon contribution. In analogy, we define the virtual 
phonon exchange contribution to the effective interaction as the contributions from ReZ) 
in Eq. (0), and the real phonon exchange contribution as that from ImD. Ignoring the 
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anisotropy factors in the phonon matrix elements for the deformation potential, it turns out 
that in the £ p ^ — > oo limit, the entire u < Qq contribution to T> 2 \ is "virtual" and the entire 
uj > c\q contribution is "real." In general, the division between the contributions are not so 
clear cut; at a given uj and q both virtual and real contributions could exist simultaneously. 

In a semiclassical transport theory, real phonon processes result in a non-equilibrium 
distribution of phonons in a coupled electron-phonon Boltzmann equation. In Appendix |A] 
we derive an expression for the real phonon exchange contribution to the drag resistance us- 
ing such a coupled Boltzmann equation approach and explicitly demonstrate its equivalence 
to the purely real phonon contribution of T> 2 \ to the drag. The appearance here of both 
virtual and real phonon contributions to the transresistivity in a single Feynman diagram is 
reminiscent of the appearance of both contributions to the quasiparticle scattering rateH~H 
in the phonon exchange contribution to the electron self-energy of a 3D electron-phonon 
system. 



B. Reciprocal space calculation 

In the following subsections we elucidate the physics of the drag for short phonon mean 
free paths. When £ ph <C £ p h, C rit, we can set the factors in the expression e to zero. The 
screened interlayer interaction is then simply 

W 21 (q,uj)~- — ~ 2 > 53 

(1 - U u (q)x(q,uj)) 2 e 2 

where we have taken the static limit of x appropriate for temperatures comparable to T BG . 
The transresistivity is then given by 

-h 2 [°° i 3 Imxi(g,Qg)Imx 2 (g,Qg) 1 r°° 2 

p2i,i^ T ^~ 2 r—] dqq —J7Z /o; =r TT du \u V 21 {q,u;)\ . 54 

8n z e z nin2k B T Jo smh (hciq/2k B T) v^e%Jo 

We have factored terms which vary slowly with respect to u out of the u integration, since 
V 2 i(q,uj) is relatively sharply peaked around uj = c\q. 

From Eqs. ( P7| ) and Q3TJ) , v^Dyi for uj w aq is given approximately by 
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v V 2 i(q,u)) ~ 

-q 2 C D p /cut(g) exp(-dq) / (k F q) , q - u/ct > 



q 2 C DP /cut (?) exp(idq - ud/2c l £ ph q)/(ik F q) ,U)jc x - q > £ 1 



(55) 

P h 



where g = |g 2 — (co>/q) 2 | 1//2 , and f cu t(<j) is a function (dependent on the form factor) which 
cuts off at q ~ L~ x . 

In the Qq > u case, the effective interaction is dominated by virtual phonon exchange, 
whereas in the u > c\q case, the integral is dominated by energy conserving intermediate 
states and the effective interaction is due to real phonon exchange. The analytic expression 
for the integral is complicated in the small region where \lo — c\q\ < Q/£ pn , and we ignore 
it in our treatment below. The smaller £ p ^ is, the wider this region becomes and hence 
the expressions derived below are not quantitatively valid for small £ p \ (below 10 4 nm for 
typical parameters in GaAs); however, the expressions seem to exhibit the correct qualitative 
behavior when compared to numerical calculations even for small £ p h_. 

Parameters of experimental systems studied to date satisfy the inequalities £ p ^ 3> d > L 
and (2/cpL) 2 ^> 1. It follows that for u near qg the integrand of the frequency integral in 
the drag resistivity expression proportional to q~ 2 ~ 2q\q — uc\\. The logarithmic divergence 
of the integral is cut off at small q by the validity limits in Eq. (|55D and at large q by the 
exponential suppression factors or cut-off functions which appear in Eq. (|55|). For the real 
phonon contribution we find 



[°° j 1 t> / m2 q^C%pQ exp(-qd/£ ph q) 
/ dco \VqD-o. (g,w) I ~ — 72 — / dq 



Iqci h, F Jq 

,3/^2 



q CdpQ 



(56) 



where 



-phj i 

r -r — 1 

ft 



' =L~\ (57) 



and Ei(x) = fl^dt exp(t)/t is the exponential integral, which has the limiting behavior 

26 



ln(x) , < x < 1; 

Ei(-x) ~ { (58) 

exp(— x)/x ,j;>1. 

The upper cut-off g max reflects the fact that it is impossible to excite phonons with a z- 
wavevector larger than L^ 1 . The q integration in Eq. (|56| ) can by physically interpreted as 
follows. The velocity component of a phonon in the ^-direction is approximately Qq/q, and 
hence the time taken to travel a distance d in the z-direction is 

Wis = qd/ciq (59) 

The exponential factor in Eq. ( |5ED is exp(— ttransQ/^ph), which is the probability that a real 
phonon emitted from one layer reaches the other. 
For the virtual phonon contribution we find 

g 3 Cg P Q f°° exp{-2dq) 



r i du\v V 2l {q^)\^ q ^^ ^ da 

JO Kp Ja v . 



F J Cin Q 

q 3 Cl P ci „ / , [2q 



- 4 Ei {-lt)- m 

where q v zmin = (g/2£ p h) 1//2 is also given by the validity limit in Eq. (55). The term in 



the exponent can be interpreted as the probability that a virtual phonon emitted by one 
layer reaches the other. The lifetime of the virtual state which is given by the energy-time 
uncertainty relation At ~ (<?/q — ^) _1 ~ q 2 /Qq. The probability is exp(— t trans / At), where 
transit time is given in Eq. (^). This gives the exponential factor and hence the distance 
dependence in Eq. (|oTi|). 

Only the sum of the two contributions will be observed in experiments. Adding the two 
terms yields 

<? 3 qCdpt,. ( qdL\ 



/ duj \is V 2 i{q,uj)\ « — 2 Ei 

JO rC rp 



g 3 QC D p v . 

F 



In (fe) , qLd/i ph « 1; 

E ph exp(-qLd/£ ph )/(qLd) , qLd/£ ph ^>l. 



Inserting Eq. fl6T|) in Eq. (p^ ) we obtain 
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Ei(-2k F Ld/£ ph ) / (f\ Imxi(g, c x q)lm X2 {q, c x q) . . 

n in2 k B T Jo q \e 4 Q ) smh 2 {hc x q/2k B T) ' 1 J 

where K\ was defined in Eq. (|44a|) and again a = 6 for A = I and a = 2 for A = t. 
The argument of the exponential integral has been replaced by the typical value £ p \j2Lkpd 
since it varies slowly over the large angle scattering region which dominates the integral 
[kp = min[fci?i, kp^}). Eq. flG2"| ) gives low T power laws for p 2 i of T 10 for A = I and T 6 for 
X = t. 

The distance dependence in the short phonon mean free path limit has now been made 
explicit; as d is increased, the transresistivity falls logarithmically when £ p \ l /2Lkpd ^ 1 and 
exponentially when £ ph /2Lk F d < 1. The dependence of the transresistivity upon tempera- 
ture, electron density and the ratio of the electron densities nxjni is given by the remaining 
integral over wavevector which essentially expresses the phase-space available for large angle 
scattering as we shall discuss in Section [VI]. 




C. Real space calculation 

In many physics problems, translational invariance and time-independence makes recip- 
rocal space calculations, like the one presented in the preceding section, simple and conve- 
nient. It is often the case, however, that the corresponding real space calculations provide 
a useful alternate language for physical interpretation. In the present case the frequency 
integral in Eq. fl54|) is equivalent to a integral over time 

/OO /"OO 
duj\u V 21 (q,uj)\ 2 = 2tt dt \u V 21 (q,t)\ 2 (63) 
-oo J — oo 

where 

voVn(q,t) ~ r dQ z F 1 {Q z )F 2 (-Q z )D(Q,t), (64) 

ZTlKp J-oo 

and the real time, retarded, phonon propagator is given by 

D(Q,t) = -2G(t) sm(uj Q t)e- t/2T ^ . (65) 
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where r p h = £ p h/ci. Since the important values of Q z are small compared to q, we can 
expand 



= c iy /q 2 + Q 2 z « Cl q + ^Q 2 Z . (66) 

iq 

Notice that the phonons have a quadratic dispersion in the Q^-direction but a linear disper- 
sion in the plane of the electron layers. The (^-integration in Eq. (|6"1D can now be performed 
if we model the density profiles as Gaussians, so that F\{Q Z )F2{— Q z ) = e~ L ^ < °&e~ %d Q z . We 
find that 

r°° c 2 a 2 G 2 r°° , ex P(~| TiTTM^) 

{^ WoVM ^^l * e ^_^_. (67) 

This expression suggests a picture in real time for the phonon mediated interaction. 
The integrand of the time integral is recognized as a wavepacket which is centered around 
one well and which broadens and decays as time evolves. When the wavepacket is broad 
enough to reach the other well, a transfer of momentum (parallel to the layers) can occur. 
The distance dependence of the interaction is determined by the time-dependence of the 
intensity of the phonon field disturbance at the other well. Since d 2 3> L 2 there will be no 
contributions to the integral in Eq. (|6~7|) until t ~ t* = qLd/ci which corresponds to the time 
it takes the wavepacket to reach the second well. For t > t* the square of the amplitude 
of the wavepacket at the second well (the integrand in Eq. fl67f)) falls off as 1/t until it is 
eventually cut off at t ~ r p h- If r p h < £*, then only the exponentially small tail of the phonon 
field impinges on the second well. Defining s = cit/qLd the time integral, Eq. ( |67D can be 
written 

2 



, . n , , l2 2Cl P Ciq 3 f 00 , 1 ( T d \ ( 

du \v Q V 21 (q, u)\ w 72 / ds -exp \-qL— s exp - 

3 kp Ji s \ £ ph j \ 



(68) 



This integral cannot be evaluated analytically, but as s > 1 the term exp(— l/2s 2 ) can 
be approximated by unity. The resulting integral then yields Eq. fl6ip . Therefore, the d- 
dependence of P21 can be interpreted as coming from the interplay of the layer separation 
and the range of the phonons which mediate the drag. 
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In the regime of short phonon mean free paths, the drag is mediated by damped phonons 
some of which spread over large distances in the z-direction. Hence, in this regime the bound- 
aries in the ^-direction will limit the effective £ p h- When boundary scattering dominates, a 
complete theory would require a realistic description of phonons scattering off the sample 
boundaries. This is in contrast to the limit £ p h ^> £ p h, C rit where the drag is dominated by a 
coupled mode which is localized on the length scale of ds and boundaries beyond this range 
are irrelevant. 



VI. NUMERICAL RESULTS 

The results presented in this section were obtained by numerical evaluation of the tran- 
sresistivity formula [Eq. (|T9|| using the complete expressions for both the effective interaction 
[Eq.(|14|)l and the screening function [Eq. (pT|)] . Although these numerical calculations are 
free of some of the approximations used in the preceding sections in order to obtain trans- 
parent analytic results, we do not expect them to be exact. In particular, corrections to the 
random phase approximation for screening must have some quantitative importance and are 
difficult to estimate reliably. 

Except where noted, the layer density was chosen to have the typical value n% = 1.5 X 
10 11 cm~ 2 which yields a Fermi temperature of Tp = Ep/kB — 60 K. The GaAs material 
parameters are taken from the literature and are given in Sec. [Ill B| .il 



A. Temperature dependence 

As mentioned in the introduction, the transresistivity increases roughly as T 2 at low 
temperatures if the effective interaction between the layers is frequency independent. De- 
viations from this behavior are indicative of retarded effective interactions. In particular, 
phonon-mediated contribution to P21 grows approximately as T 10 for the deformation po- 
tential and T 6 for piezoelectric contribution. Both cross over to a linear dependence on T at 
approximately their respective Bloch-Gruneisen temperature. In Fig. || we plot P21/T 2 as a 
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function of T. The shape of this curve is a result of an effective interaction which, for each 
q, is sharply peaked as a function of frequency around u = Qq. This peak in the frequency 
integral, which produces the dominant contribution to the phonon-mediated drag, is cut off 
exponentially by the thermal phase space factor when hciq > 2/c#T. This implies that only 
phonons with energies less than the temperature can participate in the drag. The crossover 
from T 5 or T 6 (where the piezoelectric coupling dominates) to linear T dependence of P21 
(hence the peak in P21/T 2 ) occurs when the 2kp phonons, which are responsible for most 
of the momentum transfer in the systems, can be excited. The peak is expected to occur 
near the temperature scale Tbg = KlcikpjkB which is 7.8 K for the parameters used in Fig. 
[3]. It in fact occurs at around T peak « 2.5 K ps T B g/3. For T > T BG the contribution of the 
deformation potential phonons is approximately 50 to 100 times greater than that of the 
piezoelectric phonons for the parameters used in the figure. 

We also plot the results of Eq. ([43"D and the experimental results from Ref. |j. One 



can see that the long mean free path approximation, and small layer separation, Eq. ([43]) 
underestimates the contribution of the mode. An investigation of the integrand indicates 
that this is due to a rather long non-Lorenzian tail in the uo < c\q regime which can be seen 
in the solid curve of Fig. g discussed later. 

Comparison of the magnitude of the experimental points and the theoretical curves seems 
to indicate that the effects of the collective mode have been observed. However, as we caution 



in Sec. VI D, this should not be construed as definitive proof of the existence of the mode. 



B. Density dependence 

From Eqs. (|43|) and fl62|) the transresistivity at matched densities is proportional to 
1/n 2 times an integral whose value reflects primarily the allowed phase space volume for 
exchange of phonons with large planar momentum and small perpendicular momentum. 
The maximum planar momentum transfer, q = 2kp, is proportional to the square root of 
the density. At temperatures larger than T peak , the phase-space integral should therefore 
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increase with density reflecting the larger possible momentum transfers. At temperatures 
well below T pea k, on the other hand, phonons at q = 2kp cannot be exchanged, and a higher 
density cannot be exploited. Consequently, in the case where the density is changed at a 
fixed temperature the following behavior should be observed. For low densities, the increase 
in density n permits larger momentum transfers, thus increasing the transresistivity. At 
some n, a further increase in does no good as the 2kp wavevectors cannot be excited, and 
this leads to an overall decrease in decrease in pi\. The resulting curve for p%\ has a bump 
or maximum as a function of density. The position of the bump or maximum increases with 
increasing temperature, for the reasons stated above. Fig.|] shows the density dependence 
of the transresistivity for four different temperatures, T = 1,2,3 and 4K and £ p h = lpm 
(short mean free path regime) and 1mm (coupled- mode dominated regime). 

While the general features discussed above are seen in both regimes, Fig. f| still shows a 
qualitative difference between the two. The reason is twofold. A comparison of Eqs. fl43|) 
and ( |62"D reveals that there is an extra factor of Im[x] in the integrand of Eq. (|62]) . Since 
Im[x] decreases with increasing kp, P21 falls faster with increasing density in the short mean 
free path regime than in the coupled-mode regime. More importantly, ^ p h.crit is strongly 
dependent on density: £ p h,crit oc {2k F + <2 , tf) 2 Af- Hence, increasing the density can make 
the coupled mode more and more dominant and lead to an increase of p2i- As seen in Fig. |] 
the difference in density dependence could provide a clue regarding which of the two regimes 
prevails in a particular sample. 

C. Density-ratio dependence 

The transresistivity is strongly dependent on electron-density ratio, n\jni- For tem- 
peratures well below the Fermi temperature, the maximum planar momentum transfer is 
q = 2kF,mm where kp^mm is the Fermi wavevector of the electron gas with the lowest density. 
Decreasing one density will lead to a decrease in p 2 i; a decrease also occurs if one density is 
increased because of the general \jn\Ui dependence. The resulting peak at equal densities, 
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illustrated in Fig. [|, is independent of temperature and does not occur for the Coulombic 
drag mechanism (in the absence of any plasmon effectsEUll) . This difference has been usedS 
to separate Coulomb and phonon drag mechanisms experimentally when the layer separation 
is sufficiently small that the Coulomb mechanism is of comparable importance. 

The calculations were done for both short and long mean free path regimes. However, 
the shape of the curves are so much alike that it is unlikely that a density-ratio measurement 
can be used to differentiate these two regimes in an experiment. 

D. Mean free path dependence 

In the preceding section we have explained that there is a critical £ p h beyond which a 
new collective mode appears which enhances the overall transresistivity. In Fig. |6|, we plot 
numerical results for the transresistivity as a function of the mean free path for T = 3 K and 
d = 50 nm. For £ p ^ < £ p h,crit ; we are in the regime where the electron-phonon interaction 
is separately screened in the two layers and P21 increases logarithmically with £ p h. We see 
in Fig. ^ when the mean free path exceeds £ p h, C rit and the collective mode starts to emerge, 
the transresistivity increases more rapidly before saturating at a value given roughly by 
the estimate given in Eq. ([43l) . By splitting the contributions into real (u> > ciq), virtual 
(ciq — u ^> u ) and coupled-mode (o> ~ u ), we show in the inset of Fig. | that the increase in 
P21 which occurs when £ p h exceeds ^ p h, C rit comes mainly from the coupled-mode contribution. 

In the low temperature regime where the drag experiments have been performed, £ p h 
is certainly sample dependent but should be temperature independent for a given sam- 
ple. Extrapolating the logarithmic £ p h dependence (i.e., which neglects effects of collective 
screening) in Fig. |6] for the d — 50 nm case, we find that the value of £ p h required to fit 
measured values^ of P21 (~ 12mf2) are (literally!) astronomical in magnitude. This could 
be taken as evidence for the importance of cooperative screening and collective modes in 
present experimental systems. On the other hand, our numerical calculations are dependent 
on random phase approximations estimates of Im% which are likely to require substantial 
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numerical revisions especially in the important region near q = 2k p. These corrections 
will increase Im% and could possibly boost the theoretically calculated p2i to experimental 
values without invoking the coupled electron-phonon mode. The random phase approxima- 
tion for the screening function can also lead to quantitative discrepancy. Hence, we do not 
claim that there is incontrovertible experimental evidence for the existence of the coupled 
electron-phonon mode. 

E. Layer separation dependence 

It is possible to prepare a series of double-quantum-well systems which are substantially 
identical apart from the separation between the 2D electron layers. UH Experiments on such 
a series of samples will give a more certain indication of the operative phonon exchange 
regime than experiments on a single sample. 

In Fig. [7| we plot the transresistivity as a function of the separation between the wells, 
for various phonon mean free paths. For £ p ^ small enough that the collective mode has not 
developed, P21 exhibits a purely logarithmic dependence on d, until d exceeds £ p h/(2/cj?L), 
when it begins to fall exponentially. For £ ph > ^ p h,crit, the behavior is more complex. The 
collective mode develops and p<i\ is considerably enhanced. For small d, our numerical 
results suggest that p 21 also decreases logarithmically with d, in contradiction with the d- 
independence implied in Eq. (f43"l) . Fig. || which shows the integrand of the transresistivity 
expression as a function of frequency at q = 2k p, indicates that there is a substantial non- 
Lorenzian tail for u < Qq — ujq which is not taken into account in the analysis leading to 
Eq. (f43"D. This tail contributes to the nearly logarithmic ^-dependence in the regime where 
the collective mode contributes significantly. 

Another surprising result is that in this regime, P21 does not decrease monotonically as 
d is increased. Instead there is a peak when d w y^p^crit Af- Beyond this layer separation, 
the electron-phonon collective mode of the double quantum well begins to decouple into 
two weakly coupled modes centered around each well. This is illustrated in Fig. |, where 
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we compare the frequency integrands for d < y^ p h,crit/&F and d > y ^ p h, C rit/&F- As d is 
increased further beyond ^ p h, C rit/&F the transresistivity eventually resumes its decline. 

VII. SUMMARY AND CONCLUSIONS 

Phonon exchange contributes importantly to the frictional drag resistance between 
nearby electron layers and is the dominant drag mechanism at large two-dimensional layer 
separations. In this paper we have reported on a thorough theoretical analysis of this drag 
mechanism. We find that drag includes contributions due to exchange of both virtual and 
real phonons, and that novel coupled collective modes of the 2D electron and 3D phonon 
systems can play a role depending on sample geometry and material parameters. We dis- 
tinguish two regimes based on the relationship between the phonon mean free path and a 
crossover length scale ^ p h, C rit which is typically of the order of 0.2 mm. It is possible that 
the mean free path for high quality MBE grown heterostructures can exceed £ p h, C rit at low 
temperatures. 

In the short mean free path regime the dominant drag processes at momentum trans- 
fer q have energy transfer just below c%q for virtual and just above ciq for real phonons. 
The drag rate in this case decreases logarithmically with layer separation d, until d reaches 
d a = ^ph/ZkpL ~ £ p h- The weak layer separation dependence comes nearly entirely from the 
virtual-phonon exchange contribution. For d > d a , the virtual-phonon exchange contribu- 
tion is small and the real-phonon exchange contribution decreases exponentially with layer 
separation. The real-phonon exchange contribution to the drag is consistent with expecta- 
tions based on the coupled Boltzmann equations for the electron and phonon systems, and 
the exponential fall off for d > d a can be understood in terms of the decay length for the 
disturbance of the steady state phonon system as the result of current flowing in one of the 
electron layers. 

For samples with phonon mean free paths larger than £ p h,crit a new collective mode 
involving both electronic and lattice degrees of freedom emerges below the continuum of 3D 
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phonon energies with 2D wavevectors q near Ikp. The existence of this mode enhances the 
drag. In this regime the drag also has a roughly logarithmic layer separation until d reaches 
another crossover length d& — (1 + <Ztf /2kp) /ISUfCdp- For typical samples d& ~ 0.5 /im. 
At this layer separation the collective mode separates into weakly coupled modes associated 
with the individual 2D layers. Although the the drag has a complicated and non-monotonic 
layer separation in this regime, it does ultimately decline with increasing d. These findings 
are summarized in Table (|). 

There are, unfortunately, few existing experiments on the long distance dependence of the 
phonon mediated drag resistivity. More experimental work will be needed to understand 
the range of possible behaviors and their dependence on system parameters. We note, 
however, that the preliminary experimental results by Gramila et alJH are consistent with 
a logarithmic ^-dependence. 
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APPENDIX A: BOLTZMANN EQUATION DERIVATION OF THE PHONON 

DRAG 

A sketch of the derivation of the contribution due to real emission of phonons is given 
here. This Boltzmann equation approach has be used previously to estimate the real phonon 
contribution to the drag.&S 

Assume that the distribution function of the driving layer is a drifted Fermi Dirac. Then, 
the deviation function i^i(k) = 5/o,i(k)/[(l — /o,j(k))/o,j(k)] (where f 0ti is the equilibrium 
distribution function in layer i) for a driving field eEi in the x-direction is 
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^i(k) = ^v x (k). (Al) 
The generation of nonequilibrium phonons is given by the electron-phonon coupling, and 

is 

i^W*) = ~TV ^ \9i(Q)\ 2 5 ( £ k +q - [ek + ^qD /o,i(k) (1 - / ,i(k + q) iVo(Q) x 

\ / gen k 

[Vi(k)-^(k + q) + $(Q)] (A2) 
$(Q) =«5JV(Q)/(JVo(Q)[l + iVb(0)]), (A3) 

where SN is the nonequilibrium distribution of phonons and Nq(Q) = [exp(hujQ / ksT) — 1] . 
The coupling constant #i(Q) is given by 

^(Q) = |M(Q)| 2 liMQ,)! 2 . (A4) 



Since the phonon coupling constant is small, we ignore the 3>(Q) term in Eq. (JAj) because 
it is higher order in g. 

Then, one can write the phonon generation rate (using the identity / (e)(l — fo(e+u)) = 
[/o(e)-/o(e + o;)] [N {u) + 1]), 



(^P) =~ l^(Q)| 2 iVo(g)[l + iVo(g)]^Imx(g,- Q ) 

\ / gen z B 



(A5) 



The phonon Boltzmann equation is 



d$(Q,z) xf s n] fr .., 2 qxTeE Xf , , ^ $(Q, z) ( K(C . 

v z y z = 6 ( Z ) 2 \9(Q)\ m *k B T ^ X ^ q,UJci " ' ^ ^ 

whose solution is 



2|^ 1 (Q)| 2 g :£ reE l T / 



:(Q) r ph 



(A7) 



The electron-phonon collision term in layer 2 is 

+5(^k +q - e k - fiuq) iV (Q)$(Q)}. (A8) 
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Using $(— Q) = — $(Q), the total momentum transfer to the second layer is 

fdp\ f dq qg^regi dQ z 

\dt) 2 ~ 71 J (2vr) 2 m*k B T Jo 2n 

2|gi(Q)| 2 |g 2 (Q)| 2 ex P (-d/v z (Q)r ph ) 
X MQ) 4 sinh^(^ Q /2A; B T) ^ Im * 2 ^ ^ 

Given that uq = cQ, the transresistivity can be expressed in the following form: 



(A9) 



_ E2 _ -(dp 2 , x /dt)m* 
J\ e 2 n\n 2 EiT 



-1 



8ir 2 e 2 nin 2 kBT Jo 



r 

dq q 3 / 
Jo 



00 , Imxi(g,^)imx2(g,w) 



do; 



sinh 2 (^/2A; s r) 



^(u/ctf - q^F^^/dY -~?)| 2 |M(g, ^{u/c) 2 - q 



2\\4 



exp 



In comparison, the expression for the "real" contribution to the transresistivity is 



(A10) 



P21 



~fr 2 Z -00 n 3 f°° , , rr) l2 Imxi(g,u;)Imx 2 (g,a;) 

8n z e z nin 2 k B T Jo Jo smh [huj/2k B T) 



(All) 



where 



X>2i,reai(q,w) = f°° Pi \M(q, q z )\ 2 F l (Q Z )F 2 (~Q z ) lmD(Q,u) 
J-00 Znn 



l/2r ph 



ImD(Q,cj) 

For small l/r p h the poles of the Lorenzian are 

g z = ± v / (^ + 2/(2 Q r)) 2 -g 2 



l/2r ph 



(c- QV y + g 2 ) 2 + (l/2T ph )^ 



(A12) 



±Ju 2 /c 2 — q 2 ±i- 



(A13) 



\rc 2 ^ 2 lc 2 -q 2 

As they lie close to the real axis, the Lorenzian can in general be approximated by S(uj ± 
c i\J(f + <? 2 )- The imaginary part of the pole, however, does affect the exp(— iQ z d) term 
which comes from from the difference in the phases of the form factors Fi(Q z ) and F 2 (— Q z ). 
The Q z -integration results in insertion of the imaginary part of the pole into the exponent, 
yielding 
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£>21,real(q,w) 



M q, Ju 2 /cf - q 2 



Ft ( sjutjc} -q 2 )F 2 
d 



ui 2 /cf — q 2 



exp 



cf^ 2 /cf-q 2 ' V 2QT phv /l -cfq 2 /u 2 
Substituting Eq. (gT| ) into Eq. flAllD , gives Eq. ( [MOD . 



(A14) 
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TABLES 





^ph < ^~- ^ph,crit 


^ph ^* ^ph,crit 


Physics 


Damped phonons 


Coupled e-p/i mode 


d 


d « d a 


d » d a 


d < d B 


d > d B 


Distance dependence 


Hdjd) 


^exp(-d/d a ) 


Weak (logarithmic) 


Complicated 



TABLE I. Table summarizing the distance behavior for the collective mode regime and for the 
damped phonon cases. 
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FIGURES 



FIG. 1. The Feynmann diagram corresponding to the correlation function The triangles 
(the function A) in each layer are connected by phonon propagators (wiggly lines). The dots rep- 
resent the e-ph coupling and the dashed lines are external current operators. The frequencies and 
2-dimensional wavevectors q are conserved in each vertex, whereas the perpendicular components, 
Q z and Q' z , are independently integrated over. 



FIG. 2. The two-dimensional electron density of states times the screened phonon mediated 
interaction, i^o\T>2i(2kF , oS)/ei2kpi oS)\ as a function of frequency near the longitudinal resonance 
uj ~ 2cikp. The dotted, dashed, and solid lines are for £^=0.1 mm, 0.3 mm, and 1 mm, respec- 
tively. Other parameters are d=500 A, L=200 A, and the density n = 1.5 x 10 11 cm" 2 . The 
coupled electron-phonon mode occurs when i>q Re£>2i/Re[e] < —1/2. [see Eq. (37)]. The inset 
indicates that for the given parameters the coupled electron-phonon mode develops at q = 2kp 
when £ p h si 0.5 mm. 



FIG. 3. The scaled transresistivity P12/T 2 as a function of temperature for (in increasing 
order) £ ph = 10 4 , 10 5 , 10 6 10 7 and 10 8 nm. The parameters used are n = 1.5 x 10 11 cm -2 , L = 200 A 
and d = 500 A in GaAs. The dots are data points obtained from Ref. ||, and the dotted line is 
given by Eq. (f43|). Inset: Log-log plots for the theoretical curves ^ p h = 10 4 to 10 s , showing the 
crossover from T 6 to T behavior. The dotted lines are for reference. 



FIG. 4. The transresistivity as a function of matched densities n for (a) ^ p h = 1 fim and (b) 
£ p h = 1mm. The temperatures are T = 1K,2K, 3K and 4K for dash-triple-dotted, dash-dotted, 
dashed and solid lines respectively. Other parameters as in Fig. ||. The ratio ^ p h/^ P h,crit varies in 
(b) from 0.55 for n = 0.5 x 10 11 cm~ 2 to 14.6 for n = 2.5 x 10 11 cm -2 . 
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FIG. 5. The scaled transresistivity P21T ~ 1 ni/ri2 as a function of the density ratio n\/n2 for 
(in increasing order) T = 1,2, 3, 4, 5, 6 K, n2 = 1.5 x 10 11 cm~ 2 , d = 50 nm, (a) £ p ^ = 1 /zm, and (b) 
£ph = 10 cm. Other parameters as in Fig. ^. 

FIG. 6. The transresistivity as a function of the mean free path, for d = 50, 3000 and 5 x 10 4 nm 
(solid, dotted and dashed lines, respectively). The density is 1.5 x 10 n cm~ 2 , and T = 3 K. 
Inset: Real phonon, virtual phonon and coupled-mode contributions (dotted, dashed and solid 
lines respectively) for transresistivity as a function of mean free path, for d = 50 nm, T = 3 K. 

FIG. 7. The transresistivity as a function of the well-separation for (in increasing order) 
£ ph = 10 3 ,10 4 ,10 5 ,10 6 ,10 7 and 10 8 nm, at T = 3 K, n = 1.5 x 10 n cm- 2 , L = 200 A. 

FIG. 8. The integrand as a function of to for q = 2/cp, = oo, for d = 50 nm (solid), 
d = 3 x 10 3 nm (dashed) and d = 5 x 10 4 nm (dotted). Other parameters as in Fig. 0. As the 
well-separation increases, the electron-phonon mode decouples into two weakly coupled individual 
modes, leading to the twin-peaked structure at d = 3 x 10 3 nm. At large well separations, the 
collective modes do not contribute to transresistivity, and all the momentum transfer is due to the 
real phonons. 
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